New Class of High-Energy Transients from Crashes of Supernova Ejecta with Massive 

Circumstellar Material Shells 
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A new class of core-collapse supernovae (SNe) has been discovered in recent years by opti- 
cal/infrared surveys; these SNe suggest the presence of one or more extremely dense (~ 10 5-11 cm -3 ) 
shells of circumstellar material (CSM) on 10 2 ~ 4 AU scales. We consider the collisions of the SN 
ejecta with these massive CSM shells as potential cosmic-ray (CR) accelerators. If ~ 10% of the SN 
energy goes into CRs, multi-TeV neutrinos and/or GeV-TeV gamma rays almost simultaneous with 
the optical/infrared light curves are detectable for SNe at < 20 — 30 Mpc. A new type of coordinated 
mult i- messenger search for such transients of duration ~ 1 — 10 months is required; these may give 
important clues to the physical origin of such SNe and to CR acceleration mechanisms. 



I. INTRODUCTION 

The much-anticipated era of multi-messenger astron- 
omy is coming. GeV-TeV gamma rays, which are power- 
ful tracers of cosmic rays (CRs), are detected by Fermi 
and ground-based Cherenkov detectors. Neutrinos, 
which can uniquely identify the production of hadronic 
CRs and probe dense sources, are detectable by the 
nearly-completed IceCube and the planned KM3Net 
They are useful for studying sources especially when pho- 
tons cannot escape directly. 

Violent explosions of massive stars, such as supernovae 
(SNe) and gamma-ray bursts (GRBs) may be prodigious 
neutrino/gamma-ray sources [2|. High-energy observa- 
tions should reveal their physical origin, nonthermal pro- 
cesses, and extreme environments. TeV-PeV neutrino 
detections from extragalactic sources may be possible 
for GRBs [3, 13] and certain other kinds of SNe, e.g., 
those with hidden relativistic jets [H], relativistic pulsar 
winds @ , or semi-relativistic external shocks Q • But are 
SNe with more ordinary explosions also detectable? 

In recent years, blind surveys for optical transients 
have discovered ultra-bright SNe such as SN 2008am [§| , 
2008iy ||, 2006gy Q3, 2005ap [HI and 2003ma pj, 
which are much more luminous than ordinary SNe. Their 
physical origin is not settled, but a plausible interpreta- 
tion is strong shock dissipation by collision with a massive 
(M sh ~ 1 -3OM0) circumstellar material (CSM) shell at 
R ~ 10 2 ~ 3 AU pSm, which 

is more or less analogous 
to the mechanism of type Iln SNe, though some may be 
pair- instability SNe The local rate of these ultra- 

bright SNe may be ~ 10 3-4 times larger than the local 
apparent rate of classical long GRBs [13]. But SNe with 
CSM do not have to be optically ultra-bright, and the 
existence of massive CSM shells may be even more com- 
mon. The existence of dense CSM shells is indicated from 
other SNe, e.g., SN 2006jc, 2005ip and PTF 09UJ HEIlI- 
Also, massive CSM eruptions were observed in luminous 
blue variables such as r\ Carinae [l!|, and are supported 
by the discovery of a self-obscured SN [2(| ■ 



These observations motivate us to investigate the sys- 
tem of SN ejecta crashing into massive CSM shells at 
various radii. We consider the possibility of CR acceler- 
ation at strong shocks formed by the crashes, and sug- 
gest a new class of high-energy transients that are neither 
bursts like GRBs nor persistent sources like SN remnants. 
Given their rate and timescales, new types of searches at 
intermediate time scales, with coordinated observations 
of these extragalactic SNe at < 100 Mpc, are required. 
We adopt the conventional notation Q x = Q/\G X . 



II. THE SCENARIO 

As the SN ejecta crashes into an external medium 
with density n cxt , a pair of shocks (forward and re- 
verse) are typically formed. As deceleration starts, sig- 
nificant fractions of the ejecta kinetic energy (£ e j ~ 
10 51 erg) are converted to the internal and kinetic en- 
ergy of the shocked shells. The forward shock (FS) ve- 
locity in the deceleration phase is estimated to be Vf ~ 

4000 km s" 1 ^e/5i n ex 1 /ii-^i5 3 5 2 m ^ e non-radiative limit, 
and the reverse shock (RS) velocity can be estimated 
from p C jVr ~ PoxtVj, where p is the mass density and 
V s is the shock velocity (where the index "s" indicates 
a forward or a reverse shock). In the case of ordinary 
SN shocks sweeping up the interstellar medium, deceler- 
ation becomes significant around the Sedov radius, Rd ~ 

R Scd = {§trf Z - 44 X 105 AU (M ci /M Q ) 1/3 n^/ 3 . 
If massive (M s h > M e y) CSM shells are located at 



10 5-11 



AU, with CSM shell densities of n s h 



cm , significant deceleration and shock dissi- 



pation are instead expected at Rd ~ R s h, which is 
thought to be the case for the bright SNe we consider 
here [H, [3 UH . Then, the shock velocities are evalu- 
ated via the Rankine-Hugoniot relations, and the FS will 
show a decelerating behavior self-similarly. Such violent 
eruption of massive CSM shells may be challenging for 
theories of the progenitor star, but one of the plausi- 
ble mechanisms is the pulsational pair instability, which 
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might occur at the deaths of main-sequence stars in the 
mass range ~ 95 — 130 M Q [l5l |. Ejection of smaller mass 
would be more common, and for smaller CSM mass, the 
SN shock dynamics is not as affected and the dissipated 
energy would be scaled down as oc M s \ l /{M c - ] + M s h). 

In many astrophysical shocks, sizable fractions of the 
kinetic energy seem to go into magnetic and CR energies. 
SN remnants are believed to be CR accelerators, and 
gamma rays from Galactic remnants are indeed detected 
by Fermi and air-Cherenkov detectors (e.g., 21]). It is 
expected that plasma and/or MHD instabilities can am- 
plify not only the downstream field but also the upstream 
field up to e B = 2U B /pV 2 ~ 1CT 3 - 1CT 2 (e.g., [13), 
where Ub is the magnetic energy density. If Galactic CRs 
come from SN remnants, CRs should carry e cr ~ 0.1 of 
the explosion energy, which seems consistent with obser- 
vations, 2U„/pV 2 ~ 0.1-1 HIH. In addition, CR ac- 
celeration in more dense environments (n ~ 10 8 ~ 9 cm~ 3 ) 
is suggested by detections of gamma rays from V407 
Cyg [2J] and possibly r\ Carinae [25}, which may imply 
that ~ 1 — 10% of the kinetic energy still goes to CRs. 

Analogously, we may expect such magnetic field ampli- 
fication and particle acceleration for the SN-CSM system. 
Unless the Thomson depth is tt 3> 1 — 10, one may ex- 
pect collisionless shocks, as in SN shock breakout [26[. 
The ion temperature behind the shock would become 

T? ~ 9.8 x 10 4 eV V 2 4 , (1) 

where the ion-ion collision frequency is va ~ 4 x 
10~ 4 s- 1 n w T~s /2 . Electr ons would be heated by ions 
while cooled by the Compton process and bremsstrahlung 
emission. When only the Coulomb collision is considered, 
the equilibrium electron temperature is estimated to be 

T e ~ 10 4 eV Tf 5 /5 T- 8/ 5 5 n?/ 5 (2) 

over the length of ~ u~}V s , where va > i/$ e ~ 3 x 
10~ 4 s _1 nioT e ^ 2 . Here, the photon temperature is 

T 7 S = (U°/a) 1/4 ~ 4.2 eV n^V/f in the radiation- 
dominated case. On the other hand, one may ex- 
pect that electromagnetic instabilities, e.g., i^em 
10 4 s- 1 n^ 4 , 

are faster than all relevant collision fre- 
quencies [26j ]. such as vn and fj e , so collisionless shocks 
would be formed. 

However, the shocks can be radiation-mediated when 
i? s h is sufficiently small (n is sufficiently large) and tt 3> 
c/V s (27|. For example, in the SN 2006gy-like case, we 
have 

tt = n sh a T AR sh ~ 210 n sh ,u Ai? shi i 5 . 5 , (3) 

where the radiative deceleration length for the FS, L^ cc ~ 
2-rr- ~ 1.4 x 10 15 cm nZ^^Vr} is smaller than the 

fi s h(TTV/ sn,ll /,d.o' 

width, Ai? s h- In this situation, the flow is decelerated by 
radiation in the upstream region rather than collisions or 
plasma instabilities, and CRs would not carry a signif- 
icant energy fraction, though the details are uncertain. 



Hence, we here consider CRs only when the shocks are 
not radiation-mediated. 

The CSM parameters (n s h, R s h and Ai? s h) are likely to 
have a variety of values, including cases with tt < 1 — 10, 
so that one may expect the situation where CRs are ac- 
celerated at shocks formed by the collision between the 
SN ejecta and the CSM shell. For demonstrative pur- 
poses, we suppose CSM shells of M s h ~ 1 — 30 M Q at 
R sh ~ 10 15 ~ 17 cm, which are motivated by the SN-CSM 
collision model for explaining ultra-bright SNe [3, [lj| , 
and use two representative cases: Model A (n s h = 
10 11 cm" 3 , i? sh = Ai? sh = 10 15 - 5 cm) and Model B 
{n sh = 10 7 - 5 cm" 3 , R sh = AR sh = 10 16 5 cm). As for the 
SN ejecta, for simplicity, we assume fast moving, uni- 
form ejecta with £ Ci = 10 51 erg and V^j ~ 10 km s -1 
(where one may expect that the ejecta is several solar 
masses, which is lower than the CSM mass [15]). The 
corresponding shock velocities (in the deceleration phase) 
are V f ~ 10 3 5 km s" 1 and V r ~ 10 4 km s" 1 in Model 
A, and Vf ~ lO 3 ' 7 km s" 1 and V r ~ 10 3 9 km s" 1 in 
Model B. These FS velocities are also consistent with 
observed ones (e.g., [13]). Model A is close to the case 
considered for explaining SN 2006gy [TJ, [r| , whose ra- 
diation energy and peak luminosity was £ p h ~ 10 51 erg 
and L p h ~ 10 44 erg s _1 respectively. In our setup, we 
may assume shock acceleration of CRs, except for the 
FS, since L^ cc < AR s h there. In Model B, a larger col- 
lision radius is assumed, which may be more similar to 
the case of dimmer but longer-lasting ultra-bright SNe 
such as SN 2008iy @, which had £ ph - 10 50 erg and 
Lph ~ 10 42 ' 5 erg s~ 4 , and we may consider CR accelera- 
tion at both the shocks. Note that, although we consider 
the uniform CSM shell with Ai? s h = R s h, our results are 
not much changed if the CSM shell (with Ai? s h = R s h) 
has a wind-like density distribution. 



III. PRODUCTION OF NEUTRINOS AND 
GAMMA RAYS 

CR acceleration. — When collisionless shocks are in- 
deed formed, they would accelerate charged particles to 
high energies via the Fermi acceleration mechanism, pro- 
ducing a power-law distribution, dN p /dE p oc E~ q where 
q ~ 2 (28|. The acceleration time scale is written as 

^acc = V^Z' wnere V ~ 1TW ^ n ^ ne Bohm limit. We 
assume that the magnetic field grows to a fraction of 

1/2 1/2 

the equipartition value, Bf ~ 26 G 5 n sh,n^/: 3 - 5 

(B r ~ 24 G eB / r 2 -2.5( M oj/M ) 1/2 i? 1 5 3 5 2 ^,4) I from anal- 
ogy with SN remnants. Then we have 

i{ cc ~ 2.6 x 10 4 s E p>100 TcV e-y^n^H V f ~l 5 (4) 

tl cc ~ 2.7 x 10 3 s E P}100 T C ve2%. 5 (M cj /M Q y 1/2 R% 2 5 V- 4 3 , 

which are larger than collision time scales at the thermal 
energy. 

The maximum CR energy, -E™g X , is determined by 
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comparison between the acceleration time and the cool- 
ing and dynamical time scales as well as the confinement 
condition. High-energy CRs lose their energies via adi- 
abatic losses, and pp and pj inelastic scatterings. The 
CSM density is so large that the pp reaction is efficient. 
When the pp cooling time, t pp — nK 1 
we have 



limits £™ s ax , 



771m ax 
771m ax 



81 TeV 4 / /,-2. 5 <(^/ 3 3.5 (5) 
8.6 x 10 3 TeV e)/*(M Q /M e rf /2 R% 2 5 V r \. 



l/2 T/ 3 



Here a pp ~ 10 cm is the pp cross section and n pp 



0.5 is the pp inelasticity. At large R, the adiabatic cooling 
time t a d, comparable to the dynamical time t s « R/V s , 
may become more important, 
one obtains 



When they limit ££ n s ax , 



771 m ax 
771111 ax 



3.9 x 10 4 TeV 4/ _2.5<ii^/ 3 3.5^r (6) 
1.2 x 10 5 TeV e^_ 2 jM ci /M Q ) 1/2 R- 5 3/ 5 2 V r 3 4 t^ 



■;6.5- 



At sufficiently high energies, 757 processes occur, as of- 
ten expected in extragalactic sources such as GRBs 0- 
5]. Both the photomeson cooling time (t pl ) and Bethe- 
Heitler cooling time (teii) depend on the radiation field. 
The emission from SNe is the reprocessed emission of 
ultraviolet, x-ray, and gamma-ray photons. In type Iln 
SNe, the reprocessing is expected to occur in the dense 
CSM. Calculating the spectrum requires detailed numeri- 
cal modeling including the radiative transfer [IHHII, and 
we do not treat the details of its production here. Just 
for simplicity, we adopt a black-body spectrum that is 
valid in the thermal equilibrium limit (though the ther- 
malization would be incomplete especially in Model B). 
In the case of SN 2006gy, the CSM shell is opaque even 
for Compton scattering, where stronger thermalization is 
expected and radiation should diffuse out of the shell [3] , 
and the interior temperature is set by assuming that ra- 
diation with £ p h fills a sphere with radius i? s h + Ai? s h- In 
Model B, where tt ~ 1 < c/V s , radiation can escape the 
shell in the shock crossing time. We instead estimate the 
radiation energy density to be « . ^"i* 3 , where 
L p h « £phV s I R s h- The typical energy of protons inter- 
acting with T 7 photons is Ef ~ 1.6 x 10 5 TeV T~q and 
~ 4.8 x 10 3 TeV T" 1 , respectively [6J. 

We numerically evaluated t pp , t pl and teH, and de- 
termined i?™f x by comparing £ acc to cooling time scales 
and t s Q, where we assumed £ p h = 10 51 erg in Model 
A and £ ph = 10 50 erg in Model B. We obtained £™ ax ~ 
3.2 x 10 3 TeV in Model A, and E£f ~ 5.0 x 10 3 ' TeV 
and E^ ~ 2.0 x 10 4 TeV in Model B. We found that 
both t pl and £bh are not very relevant in our cases. 

Accelerated CRs are mostly confined and produce 
mesons via inelastic pp scattering, which leads to the 
production of neutrinos and gamma rays. At i? s h, the 
efficiencies of the pp reaction during t s are estimated to 
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FIG. 1. Energy fluences of muon neutrinos from a SN crash- 
ing into dense CSM, where eb = 10~ 2 ' 5 , e cr = 0.1 and 
d = 10 Mpc are assumed. Thick and thin curves represent 
Model A and Model B, respectively (see text). The dotted- 
dashed curves show the zenith-angle-averaged ANB within a 
circle of radius 1° [3l|; we use At = 10 7 s for Model A (thick 
line) and At = 10 7 S s for Model B (thin line). 



be 0,0 



fpp ~ *f/*pp — 470 i?sh,15.5"'shaiV r 



'PP 

fcj 
J PP 



f,3.5 



(7) 



t /f ej 



13 R 



-2 

sh,15.5 



(M ej /M )V r 



r,4 



pp reaction should be efficient for typical 
< ini6 5 



Hence, the 

parameters, i? s h < 10 ie t> cm and M s h ^1 — 30 Mq, 
which are consistent with the parameters suggested for 
explaining ultra-bright SNe in the SN-CSM collision 
model [ITCH [Hj]. (In ordinary SN remnants, f pp ~ 

10 -5 (M ei /M Q ) 1/3 nl^V f ~l 5 at the Sedov time.) Al- 
though CRs may make further pions through diffusion 
even after the shock crossing, CRs would also have adia- 
batic cooling due to an expansion with ~ M C jV C j/(M e j + 
M s h) (which is the fluid velocity of the merged ejecta). 
Our estimates on the pp efficiencies are relatively conser- 
vative and adequate. 

Neutrinos. — To calculate neutrino and gamma-ray 
spectra from the pp reaction, we performed numerical 
calculations [3(| , where the CR spectrum was assumed to 
be dN p /dE p oc E~ 2 e ~( Ep/E p.° x \ normalized by the total 
CR energy, £ cr = e cl £ cj . The effects of radiative and/or 
hadronic cooling of mesons on spectra, which are relevant 
in GRBs [US] and magnetars [6[, can be neglected. 

After flavor mixing, ~ 1/6 of the pion energy is carried 
by each neutrino flavor, so the muon neutrino fluence is 

El<j> v - 6 x 10~ 2 GeV cm" 2 min[l, f pp ]e cl . -i£ C j.5idr 2 , 

.... ( 8 ) 

which agrees with our numerical results shown in Fig. 1. 
Here d — 10 Mpc d\ is the distance. The signals com- 
pete with the atmospheric neutrino background (ANB), 
so one has to know both the source direction and timing 
for detections. The SN direction can be determined well 
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by optical/infrared observations. Importantly, the high- 
energy transients from the new SN class discussed here 
have much longer duration (months-years) than GRBs 
(seconds-minutes) [U, 0] and SNe with relativistic com- 
ponents (minutes-days) [H, Q , and thus form a new type 
of neutrino sources. In the shock model for ultra-bright 
SNe such as SN 2006gy thermal photons leave the 
source when the photon diffusion time is comparable to 
the shell expansion time, where 

h _ D « i^L nshCTT „ 10? s n^V'l, (9) 

(which is consistent with the observation, £ p h ~ 10 51 erg 
and L p h ~ 10 44 ergs -1 ). For the neutrino search by 
IceCube-like detectors, we have to set a time window At, 
which is relevant to estimate the ANB. In Model A, it 
would be appropriate to use At = 10 7 s since the duration 
of the SN thermal emission is i 7 _D ~ tf ~ 10 7 s, where 
the muon yield from SN-CSM neutrinos for IceCube is 
"^,>4 TcV ~ 2. In an optically thin case like Model B, 
the SN emission time is order of 

t s ~ 5.0 x 10 7 s Rsh^V-ls (10) 

(which is consistent with £ p h ~ 10 50 erg and L p h ~ 
10 42 - 5 ergs -1 ), and we obtain N^, >2 q TcV ~ 1 for the FS 
(■Nft,>50 TcV ~ 0.2 for the RS) for this time window. For 
up-going neutrino sources, attenuation in Earth should 
be considered, but will be modest at a wide range of 
zenith angles for the most important energies [32j|. 

The rate of SNe with dense and massive CSM is un- 
certain, but a few % of all SNe may be such systems [H, 
EM H3j so that their rate within 20 Mpc is order of 
~ 0.1 yr -1 . Note that the cumulative background muon 
neutrino flux, E%$ v ~ 2.7 x 10 -9 GeV cm -2 s -1 sr -1 , 
though comparable to that from GRBs 0, H| , is less than 
the ANB up to E v ~ 300 TeV, so that we focus on de- 
tections of individual nearby explosions. 

Gamma rays. — Neutral pions lead to gamma rays that 
are interesting targets for Fermi and future Cherenkov 
telescopes such as Cherenkov Telescope Array (CTA), 
where one has to consider their interactions inside and 
outside the source. At relevant energies, there will be 
attenuation on matter (Bethe-Heitler pair-creation) and 
radiation (77 pair-creation) in the source. They are also 
attenuated by the extragalactic background light (EBL). 

In Fig. 2, the attenuated pionic gamma-ray fluxes 
are shown, taking into account gamma-ray attenua- 
tion numerically. For demonstrative purposes, the non- 
attenuated flux in Model A is also shown. Here, for sim- 
plicity, we employ ~ exp(— r 77 — tbh) for the screen 
region and ~ 1/(1 + t 77 + Tbh ) for the emission re- 
gion. For example, the Bethe-Heitler and 77 pair- 
creation depths in the CSM shell are estimated to 
be r B ^ ~ 3.2 n shj iiAi? sh: i5. 5 at ~ GeV and ~ 

3000 Tl AR shA5 . 5 at - 260 GeV T -1 , respectively. 
(Note that the photomeson and photodisintegration pro- 
cesses can also happen at sufficiently high energies. 
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FIG. 2. Energy fluxes of pionic gamma rays, corresponding 
to Fig. 1. Gamma-ray attenuation inside and outside the 
source is included; the double-dotted curve shows the intrinsic 
spectrum without attenuation in Model A. Left dot-dashed 
curves show Fermi /L AT differential sensitivities at t — 10 6 ' 5 s 
(~ U in Model A) and t = 10 7 5 s (~ t r in Model B). The 
100 hr differential sensitivity of CTA is also overlaid (right). 



When tt is sufficiently small, low-frequency synchrotron 
far-infrared emission may also increase the attenuation 
far above a TeV). Outside the source, the EBL attenua- 
tion is significant only at > 100 TeV for d ~ 10 Mpc. 

In Model A, the Bethe-Heitler and 77 attenuation 
would make it difficult to detect ~ GeV and ~ TeV 
gamma rays, respectively, although the attenuated flux 
just represents the relatively conservative case (see be- 
low). Also, though we show the differential sensitivities 
of Fermi and CTA, the integrated sensitivities over sev- 
eral energy bins are much better, which would help de- 
tection of the signal. In Model B, Tgjj <C 1 and the 
77 attenuation is negligible at < TeV, so that gamma 
rays seem detectable by Fermi for d < 20 — 30 Mpc, 
which motivates searches for ~ 0.1 — 1 yr transients via 
multi-year Fermi observations. With coordinated follow- 
up searches, ~ 0.1 — 1 TeV gamma rays may also be de- 
tected by Cherenkov telescopes such as CTA. The High 
Altitude Water Cherenkov Experiment (HAWC), with a 
larger field of view and lower sensitivity, may also be 
helpful for nearby SNe. 

The neutrino signature is quite direct and more impor- 
tant as a smoking gun of the CR acceleration, while the 
gamma-ray signature would be more complicated. While 
we are here mainly concerned with pionic gamma-ray 
emission that is the more direct hadronic signal, gamma 
rays are also produced by electrons, which come from 
muon-decay, primary acceleration, and pair-creation pro- 
cesses. When they lose energy via inverse- Compton (and 
synchrotron) processes, electromagnetic cascades can be 
induced, so that gamma-ray signals may be enhanced. 
Let us estimate the cascade effects on the gamma-ray 
emission briefly. As noted above, gamma rays of ~ 
260 GeV T - q lead to generating energetic pairs. The syn- 
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chrotron and inverse-Compton cooling times are t e _ syn 
7.7 x 10 3 s 



number of x ray photons is 



and a 



7 e IB, 2 and t„ 



2.2 x 10 2 s 77aTr 4 



respectively (while the bremsstrahlung cooling time is 
ie-brcm — 1-0 x 10 4 s tl (In 7e,3)~ )■ If interactions 
with matter could be neglected, the situation would be 
much simpler. When we consider only a one-zone black- 
body radiation field (though it may not be a good approx- 
imation), unless the synchrotron cooling becomes more 
important, we obtain an inverse-Compton cascade, which 
leads to a broken power-law energy spectrum with a high- 
energy cutoff E^ ut (e.g., [Hj|). 

In Model A, because of the large 77 pair-creation 
depth, a flat energy spectrum is expected from ~ 
0.1 GeV (£™7l0 GeV) 2 T 7i0 to E° ut ~ 10 GeV (and the 
photon index becomes ~ 1.5 at lower energies). The cas- 
caded flux is ~ 2 x 10 -8 GeV cm~ 2 s _1 (which is compa- 
rable to the intrinsic flux without attenuation), enhanc- 
ing the ~ 1 — 10 GeV flux only by a modest factor; this 
estimate was also checked via the numerical cascade cal- 
culation 3J|. However, interactions with matter cannot 
be neglected, so that even ~ GeV gamma rays are atten- 
uated mainly because of the Bethe-Heitler attenuation in 
the CSM shell and the eventual gamma-ray emission is 
expected at even lower energies, ~ 10 keV — 100 MeV, 
below the range of our interest (and note that the ra- 
diation temperature around the photosphere would be 
lower) . Then, further downgrades are unavoidable, since 
the Compton opacity becomes larger than unity for pho- 
tons, and electrons mainly lose their energies via the 
Coulomb interaction (t e _coui — L2 x 10 3 s 7e,2^Jn), 
and the bremsstrahlung process can overtake the inverse- 
Compton process. Then, the GeV-TeV emission would 
be significantly suppressed as expected above, and low- 
energy emission is largely thermalized. However, the 
CSM might be clumpy as suggested in SN 2005ip fl8j . 
where the Bethe-Heitler process might not be so signifi- 
cant depending on the unknown volume filling factor. If 
that is true, Fermi may find the signal if ~ 3% of gamma 
rays (including both the attenuated/cascaded compo- 
nents) can escape from the source, where the flat energy 
spectrum in the ~ 0.1 — 10 GeV range is expected because 
of the inverse-Compton cascade. But, since strong ther- 
malization is required in the SN 2006gy-like case, this is 
recognized as optimistic for detectability. 

On the other hand, in Model B, the 77 pair cre- 
ation is relevant only at > TeV and the Bethe-Heitler 
process is negligible, so that the attenuated flux shown 
in Fig. 2 (which is essentially the intrinsic flux with- 
out attenuation) gives the dominant contribution for 
the detectability of Fermi and CTA. Here, the inverse- 
Compton cascade is suppressed by the synchrotron cool- 
ing (i e _ syn < t e -ic) but the secondary synchrotron flux 
in the ~ 10 keV — 100 MeV range (see also the next 
section) would be enhanced by a factor. Note that the 
77 pair-creation opacity by synchrotron x rays coming 
from the generated secondary pairs is small enough. Let 
us suppose all of the pair luminosity goes into x rays 
of Ex ~ keV, though it is not true. The maximum 



47r(fl sb +Afl Bll ) 2 cBj l 

keV photon travels ~ 3.8 x 10 15 cm njjj 7 5 in a homo- 
geneous medium, so that the 77 pair-creation opacity is 

~ 6.2 X 10~° i e ,41.5-Rl _ 6.5- E; x,keV n H ) 7.5- 

Generally speaking, calculating the detailed gamma- 
ray spectrum is not easy without knowing the multi-zone 
spectral energy distribution of target photons, which de- 
pends on details of the model via the frequency de- 
pendence of the opacity and the density gradient in 
the CSM 29]. The radiation deviated from a black- 
body spectrum, including x-ray emission from thermal 
bremsstrahlung expected in the downstream region, may 
affect results. Though we defer such studies, we have 
discussed both the relatively conservative (with atten- 
uation) and optimistic (with cascade) cases, which are 
enough for our purpose of demonstrating detection po- 
tential of the gamma-ray signal. Also, one should keep 
in mind that the cascade outside the source is also ini- 
tiated, but the resulting gamma-ray emission cannot be 
detected unless the intergalactic magnetic field is weak 
enough [Hj]. 



IV. IMPLICATIONS AND DISCUSSIONS 

We have shown that SNe crashing into dense CSM 
are interesting targets for current and near-future high- 
energy detectors. Importantly, new types of coordi- 
nated multi-messenger searches are required to detect 
such 0.1 — 1 yr transients. 

CR acceleration in dense surroundings is also moti- 
vated by recent gamma-ray observations [24| [H| . Both 
detections and non-detections of SN-CSM emission will 
be useful, since physical mechanisms in such extreme en- 
vironments are uncertain, especially when ir > 1. In 
particular, neutrinos have the benefit of probing hadronic 
CR accelerators, and they can escape earlier than ther- 
mal photons, which may be delayed by diffusion. In ad- 
dition, the detections of signals would support the SN- 
CSM scenario [13; (3 rather than the pair- instability sce- 
nario [l6], useful for revealing the origin of bright SNe. 

One may expect synchrotron emission in the infrared- 
to-gamma-ray bands, as electrons of energy E e emit pho- 
tons with ~ 44 keV E 2 TcV B. But the emission should 
be reprocessed to energies < 10 keV and/or strongly 
thermalized when tt 3> 1, so we basically expect ther- 
mal emission observed from such SNe. But the syn- 
chrotron emission can be seen in the sufficiently hard x- 
ray range, if the collision happens at tt < 1, as in Model 
B. The unabsorbed energy flux is estimated to be ~ 
7xl0~ 13 erg cm -2 s _1 / syn min[l, / pp ]e cri _i£ oj , 5 i<ii~ 2 i~ 7 8 , 
where / syn < 1 is the efficiency of synchrotron cooling. 
In Model B, the synchrotron cooling is more important 
than the synchrotron self-Compton and external Comp- 
ton cooling at sufficiently high energies, so that pionic 
gamma rays are dominant in the GeV- TeV range while 
the synchrotron component, whose high-energy photon 
index is (g+2)/2 ~ 2, is relevant below ~ 100 MeV. (But 
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other cooling processes such as the Coulomb interaction 
becomes more important especially at E e < 0.1 — 1 GeV.) 
Hence, this signal will be an interesting target for the 
near-future x-ray monitor NuStar (softer x rays can be 
masked by strong thermal bremsstrahlung emission that 
seems to be observed [9(). The radio emission is sup- 
pressed by the Razin effect, free- free absorption, and syn- 
chrotron self- absorption when the collision radius is small 
enough, though it may be expected at very large radii. 

For more quantitative theoretical studies, hydrody- 
namical simulations with radiation transfer and CR back- 
reaction are desirable. But our results are enough for the 
purpose of this work, which, in part, is to motivate new 
searches starting now. The relevant quantities, £ CI , E™f^ 
and q, have uncertainties, but we could see a source up 
to d ~ 30 — 60 Mpc if e cr is larger than 0.1. The spectral 
index also affects the results. Although we adopt q = 2, 
steeper indices lead to lower muon yields from SN-CSM 
neutrinos. But harder indices may also be expected since 
the gas may be radiation-dominated and/or the shock 
may be CR-mediated [28| . Since E™f x depends on shock 
velocities, the non-uniformity of the ejecta affects the RS 
velocity and ratio of the RS dissipation to the FS dissi- 
pation [3Sj . but its pre-collision velocity distribution is 
uncertain since the ejecta may sweep the CSM before 
the collision. Also, the shock evolution may be radiative 
rather than adiabatic (36j . 

Plasma effects can modify the results, via e.g., wave 
damping by neutral particles or radiation. Especially, 
E™f* may be limited by the size of the ionization re- 
gion, since damping occurs in a time ~ — -; — — ~ 



lO 05 s n~7 5 T ~ - 4 in the neutral region 137]. Although 
the CSM gas would be initially neutral, ionization in 
the downstream and upstream region (not far from the 
shock) seems expected observationally [18[ and theoreti- 
cally [38| since the post-shock temperature is high. 

We have considered extragalactic SNe with dense and 
massive CSM shells. Such collisions may happen even 
for GRBs, where the jet breakout emission can be ex- 
pected. The high-energy neutrino and gamma-ray emis- 
sion is also possible, which would be more or less anal- 
ogous to the (sub-)photospheric emission 0,0]. In the 
Galaxy, rj Carinae is a promising candidate that showed 
violent mass eruptions [lj|. The radius of its mas- 
sive nebula is larger than our typical values that are 
required for bright SNe. But the neutrino detection 
seems possible if the star explodes, since iV^ >i%v ~ 
2 x lO 5 f cr , 5 o(M sh /lOAf )i?^ 2 17 V;73 1 5 (where f„ < 1). 
For smaller CSM mass, though the radiation might push 
the CSM shell, SN dynamics are not largely affected as 
in ordinary SNe [39|, where £ cr , E™£* and are much 
smaller. Detections would be challenging, but it may be 
interesting for a Galactic event. 

After this work was submitted and put onto the arXiv 
(|arXiv:1012.2834| , we became aware of Ref. 03], which 
is closely related to ours, and which supports our claims 
that these unusual SNe are interesting and that their 
high-energy emission is an important probe. We thank 
K. Ioka, C. Kochanek, C. Rott and R. Yamazaki for dis- 
cussions. This work is supported by JSPS and CCAPP 
(KM), E. C. Howald Presidental Fellowship (BCL), Sloan 
Fellowship and NSF Grant AST-0908816 (TAT), and 
NSF CAREER Grant PHY-0547102 (JFB). 
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